Superlattice with hot electron injection: an approach to a Bloch oscillator 
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A semiconductor superlattice with hot electron injection into the miniband is considered. 
The injection changes the stationary distribution function and results in a qualitative change 
of the frequency behaviour of the differential conductivity. In the regime with Bloch oscil- 
lating electrons and injection into the upper part of the miniband the region of negative 
differential conductivity is shifted from low frequencies to higher frequencies. We find that 
the dc differential conductivity can be made positive and thus the domain instability can be 
suppressed. At the same time the high-freguency differential conductivity is negative above 
the Bloch frequency. This opens a new way to make a Bloch oscillator operating at THz 
frequencies. 

Miniband electron transport in semiconductor superlattices is an important subject of experi- 
mental and theoretical investigations beginning from the pioneering paper of Esaki and Tsu [|l| . If a 
static electric field E s is applied along the superlattice axes, electrons begin to move in accordance 
with the semiclassical Newton's law (neglecting scattering) 
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where p is the momentum along the superlattice axes. If the energy gap A^ between the first and 
the second miniband is large enough, Aq S> edE s (d is the superlattice period) and the scattering 
rate v is small, hv < edE s , then electrons oscillate inside the first miniband with the so-called Bloch 
frequency Q = edE s /Ti. The quasiparticle energy e p = (A/2)(l — cos(dp/H)) and the quasiparticle 
group velocity along the superlattice axes v p = de p /dp = (Ad/2h) svn{dp/K) are periodic functions 
of time (A is the miniband width). We assume that the condition of semiclassical approximation 
A » edE s is fulfilled. 

In semiconductor superlattices the Bloch frequency can be varied by the electric field up to 
1-10 THz, so it is a great challenge to make a tunable Bloch oscillator in the THz frequency range. 
However, in most cases the phases of Bloch oscillating electrons are random and high-frequency 
emission is noise-like. The amplification of an external high-frequency signal is possible due to a 



2 



phase bunching of electrons. In 1972 S.A. Ktitorov, G.S. Simin, and V.Ya. Sindalovskii |2j showed 
that the linear differential conductivity cte(uj) of a biased superlattice is negative up to the Bloch 
frequency (see Fig. 0(a), broken line). Later A. A. Ignatov, K.F. Renk, and E.P. Dodin |3| calculated 
the finite-field ac response of a biased superlattice and showed that the amplification efficiency, i.e. 
the ratio of power absorbed by a THz field and the power delivered from a bias source, is larger 
at larger ac fields. That is what one needs to make an oscillator work at u ~ 0. However, o~e{w) 
is negative also at u — > 0, and this dc negative differential conductivity (NDC) leads to a low- 
frequency instability of space-charge waves and to the formation of domains |H| . This effect can 
be used itself for a generation of microwaves by moving domains with frequencies above 100 GHz 
||, but the efficiency of a THz radiation at Bloch frequencies cannot be large in the inhomogeneous 
regime. For a Bloch oscillator to operate the domain instability should be suppressed. It was shown 
by Yu.A. Romanov, V.P. Bovin, and L.K. Orlov j?J and by H. Kroemer that in the nonlinear 
regime with high enough high-frequency current the dc differential conductivity is positive, while 
the large-signal high-frequency differential conductivity still remains to be negative, consequently, 
the steady-state operation of a Bloch oscillator in large-signal regime can be achieved. However, 
the question of "device turn-on" is not answered. Note also, that the nonlinear regime can be 
unstable ||. 

We present here a new approach to this long-standing problem. The key idea is to inject 
electrons into the upper part of the first miniband (note the discussion of other new effects in a 
superlattice with hot electrons in Jl(|). We calculate the nonequilibrium stationary distribution 
function, and both linear and nonlinear (large ac current) differential conductivities with the help 
of a one-dimensional kinetic equation and 3D Monte-Carlo simulations. We show that the low- 
frequency domain instability is suppressed and the region of negative differential conductivity is 
shifted from low frequencies to higher frequencies. Thus the self excitation of high-frequency 
oscillations becomes possible. At a finite level of high-frequency current this effect is supported by 
dynamic domain suppression found by Romanov and Kroemer. 

We restrict our analytical consideration to the one-dimensional problem, where the dynamics 
of the system can be described by a semiclassical Boltzmann equation for the distribution function 
f(p,x,t). The distribution function determines the particle density 

/■nh/d 
f(x,p,t)dp, (2) 
-nh/d 

and the current density 

/izh/d 
v p f(x,p,t)dp. (3) 
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In a one-miniband approximation with a Bhatnagar-Gross-Krook (BGK) scattering integral, 
which conserves the number of particles and makes it possible to describe space-charge effects [11], 
we have 

%+->g+<->{'-&) +a <* (4) 

where 

n d ( A dp\ 

h = 2nhI (A/2T) eXP l^ COS Tj (5) 

is the equilibrium distribution function of a non-degenerate electron gas and no is the equilibrium 
particle density. S(p) is the source of hot particles. For simplicity (and to conserve the total 
charge) we choose a "narrow" source of the simplest form 

S(p) = Qd(p-p')-^-f s (p), (6) 
n 

where f s {p) is the stationary (static and homogeneous) solution of (||). Q is the injection rate of 
hot electrons and p' is their momentum. 

In the homogeneous case (constant electric field E s ) and S(p) = the solution of Eq. (||) is 

27rhl (Aj V \2TJ v + im 
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where 99 = dp/Ti is the dimensionless momentum, /; is the modified Bessel function of l-th. order. 
Note, that due to the one-miniband approximation f(x,<p,t) is a periodic function of ip with a 
period of 2ir. 

In presence of the hot electron source (||) the stationary homogeneous distribution function 
f s ((p) can be represented as f s = fE J r f and f'((p) can be obtained from the equation 

n— = - zH /+_ 6{tp-ip) JEW- (8) 

u(f \ no / n no 



Thus we find 



n d^ Qe^ 
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2ir% i v + "H^ + 
where we have introduced the parameter of nonequilibrium 

Q 



I + 



(9) 



(10) 



Distribution functions f s {(fi) at different injection rates are shown in Fig.^(a) for a superlattice 
with A = 16 meV, d = 5.06 nm, v = 0.5 • 10 13 sec" 1 , E a = 18 kV/cm, O = 1.38 • 10 13 
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FIG. 1: Stationary distribution functions in the miniband: (a) analytical results for fo(<p) (dots); /e(<^) 
(broken line); f s (<f) at 77 = 0.036 (thin line) and 77 = 0.072 (thick line); (b) 3D Monte-Carlo simulation 
without injection (crosses); at r\ = 0.036 (triangles) and i] = 0.072 (circles); is/fl = 0.36, A/2T = 0.3, 
<p' = 0.9w. 



sec -1 , at room temperature. We also calculated the distribution function by using a one-particle 
Monte-Carlo method Jl^], where the motion of an electron described by the Eq. (|l|) is interrupted 
by three-dimensional scattering events (on acoustical and optical phonons), randomly distributed 
in time with the average frequency v. The averaged frequency for each scattering process was 
calculated from quantum mechanical perturbation theory taking into account the superlattice 
energy dependence and the material parameters of bulk GaAs. The injection is modeled in the 
same manner by the random change of the momentum to p' with the averaged frequency Q/uq. 
The results of our Monte-Carlo simulations are presented in Fig.[I](b). The comparison with the 
analytical results shows that in this parameter range the BGK scattering integral in the one- 
dimensional Boltzmann equation (Q) adequately models the scattering in the 3D system. 

The next step is to find the high-frequency conductivity in this nonequilibrium state. We 
consider perturbations with frequency u> and wave- vector k of the form E(x, t) = E s +E u) ^e~ tujt+tkx , 
f(x,tp,t) = f s (ip) + fuj.k{ L P)&~ %ut+tkx and solve the linearized kinetic equation 
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FIG. 2: Differential conductivities (a) analytical results: <te(oj) (broken line); cr(u>) at rj = 0.036 (thin line) 
and r\ = 0.072 (thick line); (b) 3D Monte-Carlo simulation: <r E (ui) (crosses); a nl {uj) at r) = 0.036 (triangles) 
and r? = 0.072 (circles); v/n = 0.36, A/2T = 0.3, if' = 0.9tt. 
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-(u + iv)/Q, k 



Ad 
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For the differential conductivity we find (following the method developed in Jll]] ) 



ct(uj, k) 
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where we introduced plasma frequency uj p = y^^~i e> ls the interlayer dielectric constant, Eq 
is the vacuum dielectric constant. This is the general expression for any nonequilibrium function 
f s (ip)- With the injection source (|6|) and f s {f>) given by @ we obtain 



(J E (UJ,K) 



a(uj, k) = cte{u, k) + a'(u, k), 

oo / 
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<te(u;, k) was obtained previously by Ignatov and Shashkin [11]. 

Now let us consider the conductivity in long- wavelength limit <t{uj) = a(ui,K — > 0). Taking into 
account that J n (x) — > for x — > and n > 0, we obtain 

Jl i/ 2 (n 2 -t/ 2 + «/u;) 

(Jfi^J = <7q — - — — — (15) 

i (f + ii^)[(o; + ivy — \V\ 
which was first obtained by Ktitorov et al., and 



x < (Q — vv + iviS) cos if' — Q,(u + iv + iz>) sin <y?' 



(£ 2 + ft 2 )[(oj + ^) 2 -f} 2 ] 

7i z^ 2 (P + v) + ^(z>^ - Sl 2 )(iw - i/) 



(i/ 2 + n 2 ) 



(16) 

where o"o = , v = v + nfi, the arguments A/2T of Bessel functions are omitted. 

v 

The differential conductivity with hot electron injection is shown in Fig. ^(a) (solid lines) to- 
gether with the conductivity without injection (broken line). At strong enough injection (thick 
solid line) the hot differential conductivity is significantly different from the cold one: now NDC 
takes place above the Bloch frequency, and in the entire low-frequency part it is positive. The 
results of the corresponding Monte-Carlo simulations at large ac signal (E^ = 12 kV/cm) are 
shown in Fig.|2|(b). They confirm our main conclusion about the "inversion" of the differential 
conductivity in a nonequilibrium state and demonstrate that at high fields the effect can be even 
stronger. 

Note, that in the regime with Bloch oscillations (fi > v) the important parameter r] is deter- 
mined by the competition between the pumping frequency vt = Q/uq and Bloch frequency 0. The 
desirable effect takes place when r] is larger then some critical value r\ c . Large r\ can be achieved at 
low temperatures with a low equilibrium density uq. On the other hand, at high (room) temper- 
atures, when T > A, r\ c becomes smaller, because the equilibrium distribution function is flat in 
the miniband and therefore the effect of pumping is more pronounced. That gives hope to observe 
the effect at low as well as at room temperatures. 

An important question is the choice of the injection momentum p'. Our calculations show that 
the effect takes place if particles are injected at high enough energies, e p > A/2. The injection with 
p' and —p' leads qualitatively to the same result, but at positive p' the effect is stronger and exists 
in a larger frequency range. It can be easily seen from Eq. @ that W(u -> 0) oc — cos ip' (at 
Q 2 > v(y + rjfl)) and $ta'(u> ^> $7) oc — sin ip', such that the low-frequency instability is suppressed 
equally by p' and —p' injection, but the high-frequency instability with u> ^> exists only at <p' > 0. 

In conclusion, we have shown that hot particle injection into the upper part of the miniband 
leads to a qualitative change of the frequency behaviour of the differential conductivity. The 



7 



region of NDC is shifted from low frequencies (below the Bloch frequency) to higher frequencies. 
Consequently, the domain instability is suppressed. To further support this result we considered 
the instability of space-charge waves with finite wave vectors, the spectrum u(k) of which is given 
by the solution of the dispersion equation 

e^,k)=e' + i^^ = 0. (17) 

EqLO 

In the high- field limit (Q ^> u, u p ) we found that a low- frequency drift-diffusion mode with u(k) <C 
f2, which is unstable without injection, now becomes damped. This means, that the homogeneous 
state becomes stable, and a superlattice in that state can be used as high-frequency oscillator if 
we couple it to some external circuit (resonator or wave-guide). 

The narrow hot particle source (^) may be realized by tunneling from a metal or doped semi- 
conductor fl3| , [l4| ], or by tunneling from a narrow miniband superlattice. In the case of tunneling 
from the end of a superlattice "from layer to layer" our approach is valid for an active region with 
the length of the order of the Bloch length Ib = A/2eE s S> d. The consideration of an inhomoge- 
neous distribution function in a real superlattice with leads requires more sophisticated approaches 
|i~5| , filf l, which link Bloch oscillator and quantum cascade laser physics. 

Qualitatively the same results are obtained for the optical excitation when S{p) oc [5(p — p') + 
5(p + p')]. The consideration of optical pumping is the subject of a future publication. We think 
that the prospective way to the practical realization of the considered effect is the combination 
of injection (or optical excitation) of particles into the second miniband with relaxation from the 
second miniband into the upper part of the first active miniband. 
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